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Abstract 

Deep-inelastic ep scattering data taken with the HI detector at HERA and corresponding 
to an integrated luminosity of 106 pb“^ are used to study the differential distributions of 
event shape variables. These include thrust, jet broadening, jet mass and the C-parameter. 
The four-momentum transfer Q is taken to be the relevant energy scale and ranges between 
14 GeV and 200 GeV. The event shape distributions are compared with perturbative QCD 
predictions, which include resummed contributions and analytical power law corrections, 
the latter accounting for non-perturbative hadronisation effects. The data clearly exhibit the 
running of the strong coupling as{Q) and are consistent with a universal power correction 
parameter ao for all event shape variables. A combined QCD fit using all event shape 
variables yields as{mz) = 0.1198 ± 0.0013 toS and oq = 0.476 ± 0.008 
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1 Introduction 


This paper presents a study of the hadronie final state in deep-inelastie scattering (DIS), ep —> 
eX, using event shape variables. Event shapes probe the energy-momentum flow and are sensi¬ 
tive to perturbative QCD (the hard scattering) and non-perturbative QCD (hadronisation). Com¬ 
pared to some other hadronic final state observables, such as jet cross sections, which employ 
a limited part of the phase space only, the event shape variables include the full statistics and, 
by definition, are rather insensitive to hadronic energy scale uncertainties. With DIS at HERA 
a wide range of the scale Q is available for analysis in a single experiment, where Q is the 
four-momentum transferred by the exchanged boson from the electron to the proton. 

Studies of event shapes require an appropriate treatment of hadronisation effects, which 
often are corrected for using phenomenological models built into event generators. Alterna¬ 
tively, these effects have been treated by power corrections 0{1/Q) as calculated analytically 
from first principles by extending perturbative methods into the non-perturbative regime [1-3]. 
Within this framework, the event shapes are described by the strong coupling constant as and 
an effective coupling parameter ao for the hadronisation corrections. Measurements of event 
shapes thus represent a sensitive test of the power correction approach, and they allow ag to be 
determined. 

Previous analyses of mean values of event shape variables in DIS, as performed by the 
HI [4] and ZEUS [5] collaborations at HERA, support the power correction approach. An 
observed large spread of values obtained for as, however, indicates that higher-order QCD 
corrections cannot be neglected. Support for power corrections as an appropriate description of 
the hadronisation is also provided by analyses of mean values and distributions of event shape 
variables measured in experiments [6-8] . 

Early studies of event shape distributions in deep inelastic scattering revealed that the fixed- 
order QCD calculations available at the time were insufficient for describing the data [9]. 
Higher-order corrections have since become available in the form of soft gluon resummed cal¬ 
culations matched to NEO matrix elements [10]. This puts the study of event shape distributions 
and of the interplay between perturbative and non-perturbative QCD in the description of the 
hadronic final state on a new quantitative level. 

In this paper the distributions of five event shape variables are studied. A QCD analysis 
is performed based on the resummed and matched calculations supplemented by power cor¬ 
rections. This leads to determinations of the power correction coefficient ao and of and its 
dependence on the scale Q. 

2 Event Shape Variables 

The aim of this analysis is to study event shape variables within QCD as a function of the 
relevant hard scale, which in DIS is taken to be the four-momentum Q of the exchanged boson. 
In addition to the hadronic flow from the hard scattering there are also hadrons in the final state 
stemming from the proton dissociation, which occurs at much lower scales of about the proton 
mass. Therefore it is necessary to separate this proton remnant from the hard scattering part of 
the event. 
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2.1 The Breit Frame 


In the quark parton model the separation of the proton remnant from the hard seattering is 
elearest in the Breit frame of referenee, defined by 2xp + g = 0, where x is the Bjorken sealing 
variable, p the momentum of the proton and q the momentum of the exehanged boson. The 
z axis in the Breit frame is defined to eoineide with the proton-boson axis, the proton moving 
in the +z direetion. Partieles from the remnant are almost eollinear to the proton direetion, 
henee the hemisphere defined by pseudorapidity* p > 0 is labelled the remnant hemisphere. 
In eontrast, in the quark parton model the struek quark populates only the eurrent hemisphere 
(V < 0). 

Higher order proeesses generate transverse momenta in the final state and may even lead to 
partieles from the hard subproeess leaking into the remnant hemisphere. Still, without know¬ 
ledge of the detailed strueture of the hadronie final state, the Breit frame allows for optimal 
separation of the eurrent region from the proton remnant. All event shapes referred to in this 
paper are defined using partieles in the eurrent hemisphere only. 

The kinematie quantities needed to perform the Breit frame transformation are ealeulated 
using the eleetron-sigma method [11,12]. The virtuality of the exehanged boson is reeon- 
strueted using the energy and polar angle of the seattered eleetron, and the inelastieity y is 
determined employing in addition the energy and longitudinal momentum of all hadronie ob- 
jeets measured in the laboratory frame. This method results in good experimental resolution 
and is relatively insensitive to initial state QED bremsstrahlung. 


2.2 Definition of Event Shape Variables 

The event shape variables studied in this paper are those for whieh the ealeulations of power 
eorreetions are available. They are defined as follows. 

The Thrust variable r measures the longitudinal momentum eomponents projeeted onto the 
boson axis. It is defined as 

r = 1 - T with T = ^ . (1) 

Z2h \Ph\ 

The variable tc ealeulates the Thrust with respeet to the direetion ut whieh maximises the 
sum of the longitudinal momenta of all partieles in the eurrent hemisphere along this axis. It is 
defined as 


Tc = i- — Tc with Tc 


max 

fiqp 


Y.h \vh • Tt 
T.h\Ph\ 


( 2 ) 


The Jet Broadening B measures the sealar sum of transverse momenta with respeet to the boson 
axis 


_ '^h \Pt,h\ 

_ 2Ejpi.r 

'The pseudorapidity is defined as 77 = — In tan(0/2) with 6 the polar angle with respect to the z axis. 


(3) 
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The squared Jet Mass p is normalised to four times the squared sealar momentum sum in the 
eurrent hemisphere 


P = 


(2EJpll)^ 


(4) 


In the following the symbol po is used, whieh indieates that in the above definition the hadrons 
are treated as massless, replaeing the energy by the modulus of the 3-momentum \ph\. 


The C-Parameter is defined as 


C = 


^J2h,h' IPhWPh'lcos"^ Ohh' 

2 ~~ 


where 9hh' is the angle between partieles h and h'. 


(5) 


In Eq. 1-5 the momenta are defined in the Breit frame and the sums extend over all partieles 
in the eurrent hemisphere. 

An event is only aeeepted if the energy in the eurrent hemisphere exeeeds some value eum- 
This is neeessary to ensure the infrared and eollinear safety of the observables, beeause higher 
order proeesses may lead to event eonfigurations in whieh the partons are seattered into the 
remnant hemisphere and the eurrent hemisphere is eompletely empty, exeept for arbitrarily soft 
emissions. In the analysis the events are required to fulfill the eondition 




^ *^lim — Q/10 , 


( 6 ) 


as part of the event shape definitions. The preeise value of this eut-off turns out not to be 
erueial [4]. 

The event shape variables may be distinguished aeeording to the event axis used. The def¬ 
initions of T and B employ momentum veetors projeeted onto the boson direetion, while the 
others do not, like their eounterparts in e+e“ reaetions. Explieit use of the boson direetion im¬ 
plies sensitivity to radiation into the remnant hemisphere through reeoil effeets on the eurrent 
quark [10]. 

Throughout the paper the symbol F will be used as a generie name for any event shape 
variable. Note that for all variables F tends to zero in the ease of quark parton model reaetions 
negleeting hadronisation effeets (small values of F eorrespond to peneil like eonfigurations of 
the hadronie final state). Theoretieal ealeulations of event shape distributions and means are 
diseussed in seetion 5. 


3 Experimental Technique 

The data were eolleeted with the HI deteetor at HERA during the years 1995 - 2000 and 
eorrespond to an integrated luminosity of Eim = 106 pb“^. The eollider was operated with 
eleetrons or positrons of E^ = 27.6 GeV and protons of Ep = 820 GeV or Ep = 920 GeV, 
yielding eentre-of-mass energies ^/s of 301 GeV and 319 GeV, respeetively. Eor the present 
study, three data samples are used: 


6 




• e+p, Vs - 301 GeV, = 30 pb"^ (1995-1997); 

• e~p, ^/s ~ 319 GeV, £int = 14 pb“^ (1998-1999); 

• e+p, Vs - 319 GeV, Am = 62 pb'^ (1999-2000). 

The identification of neutral current DIS events is based on the reconstruction of an event 
vertex and of the scattered electron^ in the central tracker and the calorimeter. 


3.1 HI Detector 

A detailed description of the HI detector can be found in [13,14]. The most important detector 
components for the present analysis are the liquid argon (LAr) calorimeter and the central track¬ 
ing system. HI uses a right-handed coordinate system with the axis along the beam direction, 
the -\-z or “forward” direction being that of the outgoing proton beam. The polar angle 6 is 
defined with respect to the 2 ; axis. 

The LAr sampling calorimeter (4° < 6* < 154°) consists of lead/liquid argon electromag¬ 
netic sections and stainless steel/liquid argon sections for the measurement of hadronic energy. 
An in situ calibration provides energy scales. The electron energy scale uncertainty in the LAr 
calorimeter varies between 0.7% and 3% [15]. The hadronic energy measurement is performed 
by applying a weighting technique to the electromagnetic and the hadronic components of the 
energy deposition in order to account for the unequal response to electrons and hadrons. The 
systematic uncertainty on the hadronic energy scale amounts to 2%. In the backward region 
(153° < 0 < 177°) energy is detected by a lead/scintillating fibre Spaghetti-type calorime¬ 
ter [16]. 

The central tracking system (25° < 9 < 155°) is located inside the LAr calorimeter and 
consists of drift and proportional chambers. The chambers and calorimeters are surrounded 
by a superconducting solenoid providing a uniform field of 1.15 T inside the tracking volume. 
The scattered electron is identified by associating tracking information with the corresponding 
electromagnetic cluster in the LAr calorimeter. The electron scattering angle is known within 
3 mrad. 

For the present analysis the hadronic final state is reconstructed from combined objects, 
built from calorimeter clusters and tracks, using an energy flow algorithm which ensures that no 
double counting of energy occurs. Compared to clusters alone the combined objects improve 
the reconstruction of low momentum particles. 


3.2 Event Selection 

Several quality cuts are applied to the data. A scattered electron has to be found in the LAr 
calorimeter with a reconstructed energy 77' exceeding 11 GeV, which ensures a trigger effi¬ 
ciency above 98%. The z position of the event vertex has to be reconstructed within ±35 cm 

^Here and in the following “electron” is used to refer to both electron and positron unless explicitly stated 
otherwise. 
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of the nominal position of the interaetion point, whieh reduces contributions from beam in¬ 
duced background. Non-ep background is further reduced by requiring an event timing which 
matches the HERA bunch crossing. To suppress badly measured events the missing trans¬ 
verse momentum has to be below 15 GeV. The total longitudinal energy balance must satisfy 
40 GeV < E!i{l — cos 9i) < 70 GeV, where the sum runs over all detected particles. This 
reduces photoproduction background and initial state photon radiation. Such QED radiative 
effects, backgrounds and poorly reconstructed events are further suppressed by demanding that 
the calorimetric energy measurement of the lepton be consistent within 10% with that derived 
from the double angle method [17,18]. 

The kinematic region covered by the analysis is defined by ranges of boson virtuality 
and inelasticity y: 

196 <Q^ < 40,000 GeV^ , 

0.1 < 2/ < 0.7 , 

which are reconstructed using the electron-sigma method [11,12]. 

In total about 108, 000 events satisfy the selection criteria. The PYTHIA [19] program is 
used to estimate the background from photoproduction events, in which a hadron in the EAr 
calorimeter is misidentified as an electron. This background is found to be negligible in all bins 
of the various event shape distributions. 


3.3 Correction to the Hadron Level 

The correction of the data for limited detector acceptance uses the simulation program RAP- 
GAP 2.8 [20] with parton showers and string fragmentation as implemented in JETSET [21]. 
The parton density functions (pdfs) of the proton are taken from the CTEQ5E [22] set. The 
data are corrected for QED radiation effects using the HERACEES [23] program. Bin-to-bin 
correction factors are determined from the Monte Carlo event samples passed through a detailed 
simulation of the HI detector and subjected to the same reconstruction and analysis chain as the 
data. 

The effects of limited detector resolution are corrected for in a separate step. In the bins used 
for the QCD analysis presented in section 5 the purities are typically 30—50%. To correct for the 
corresponding bin correlations the covariance matrix is determined with an iterative Bayesian 
unfolding method [24]. The data are unfolded to the level of hadrons in order to compare with 
resummed calculations supplemented by power corrections for hadronisation. The correction 
procedure is performed separately for the three data samples. 


3.4 Experimental Uncertainties 

Several studies are carried out to estimate the experimental systematic error by using alterna¬ 
tive settings or assumptions in the Monte Carlo programs. The correlations of these changes 
between bins of an event shape variable are taken into account unless stated otherwise. 
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The event kinematies and therefore the boost to the Breit frame of reference depends strongly 
on the momentum of the reconstructed electron. Thus the electromagnetic energy scale is varied 
by its uncertainty of ±(0.7% — 3%), depending on the position of the electron cluster within 
the LAr calorimeter. In addition, the polar and azimuthal angles of the electron are changed by 
±3 mrad each. 

The event shape variables are by definition insensitive to variations of the overall hadronic 
energy scale. The effects of the uncertainty on the hadronic intercalibration are investigated 
by shifting the calibration constants of neighbouring calorimeter regions with respect to each 
other [25]. Since the final event shape distributions average over all calorimeter regions, the 
resulting uncertainty on the event shape variables is small. 

The model dependence of the correction procedure is estimated by replacing RAPGAP 
with the DJANGOH 1.2 [26] event generator, which employs the colour dipole model of ARI¬ 
ADNE [27] to simulate higher order QCD radiation. The model uncertainty is estimated as the 
difference between the results obtained with the two Monte Carlo samples. One half of this 
is treated as uncorrelated between the bins. The other half is assumed to be fully correlated 
between bins. 

An estimate of the possible intrinsic bias from the unfolding procedure is obtained by un¬ 
folding Monte Carlo event samples with themselves. The residuum between the correct hadron 
level and the unfolded result is taken as the unfolding error, which is assumed not to be corre¬ 
lated between the bins of the distributions. 

All systematic errors from different sources are added in quadrature, the model uncertainty 
and the electromagnetic energy scale uncertainty being the largest individual contributions. 


3.5 Combination of Data Sets 

The combination of separately unfolded data sets proceeds in two steps, first combining the 
positron data from the two different centre-of-mass energies and second combining the positron 
and electron data. The distributions of the two positron data sets are compatible with each other 
within errors. They are combined by calculating the luminosity weighted averages for all bins 
of the distributions. 

In general, the event shape distributions of the and e~ data are in very good agreement 
with each other, though some discrepancies are observed at the highest Q scales. Differences 
in the event shape distributions between the e+ and e~ data sets are expected because of Z 
exchange contributions to the cross section [28]. Unfortunately, only the 7 exchange component 
is accounted for in the calculations of event shape distributions which are fitted to the data. The 
e~^ and e~ event shape distributions are thus averaged, weighted with the corresponding cross 
section such that most of the Z contribution cancels. Due to the smaller integrated luminosity, 
the e~ data increase the statistical error of the final spectra. However, this effect is partly 
compensated by the larger cross section for e~p scattering at high Q. 

The resulting mean values of Q and x are slightly modified by the combination procedure, 
giving the values listed in Table 1. The average centre-of-mass energy of the combined set is 
a/s = 316 GeV. 
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# of Q bin 

1 

2 

3 

4 

5 

6 

7 

Q Interval/GeV 

[14,16] 

[16,20] 

[20,30] 

[30,50] 

[50,70] 

[70,100] 

[100,200] 

(Q)/GeV 

14.9 

17.7 

23.8 

36.9 

57.6 

80.6 

115.6 

(x) 

0.00841 

0.0118 

0.0209 

0.0491 

0.116 

0.199 

0.323 


Table 1: Q intervals, mean Q and mean Bjorken x for the seven Q bins of the analysis. 

4 Event Shape Measurements 

The normalised event shape distributions at the hadron level for thrust, jet broadening, jet mass 
and the C-parameter are shown in Figs. 1 and 2 over a wide range of {Q) = 15 — 116 GeV. 
The data points represent integrals over the bins not applying bin eentre correetions. Except for 
the highest Q bins, the precision of the measurements is not statistically limited. 

For each variable the shape of the spectra changes considerably with increasing Q, becom¬ 
ing narrower and evolving towards low values. The strong Q dependence of these spectra is 
characteristic of QCD. The results of the fits described in detail below are shown for compari¬ 
son. 

The present paper focuses on differential distributions of the event shape variables. How¬ 
ever, in order to allow for a comparison to previous analyses the mean values of these variables 
are also determined. The results are shown in Fig. 3 as a function of Q. A steady decrease of 
the means with rising Q is observed. The results of the present analysis are in agreement with 
those previously obtained [4]. 

5 QCD Analysis 

5.1 Phenomenology 

The QCD calculations used here contain a perturbative part (pQCD) dealing with partons and 
in addition use power corrections (PC) to describe the hadronisation. The perturbative part 
is made up of two contributions: fixed order terms calculated to next-to-leading order (NFO) 
in the strong coupling constant and resummed terms in the next-to-leading-logarithmic (NFF) 
approximation. 

The NFO contribution consists of the first two terms of the perturbative expansion in 
and, for any event shape variable F, has the form 

1 

-= CiiF, Q)asifir) + C2(F, Q)al{fir), (V) 

cr Clx^ 

where jir is the renormalisation scale, chosen to be Q. In order for the truncated series to be a 
good approximation to the exact solution, as{Q) needs to be small, i.e. Q should be large. 

The coefficients ci and C 2 can be calculated from the matrix elements of the hard scattering 
and the parton density functions of the proton. The fixed order coefficients in the MS scheme are 
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Figure 1: Normalised event shape distributions eorreeted to the hadron level for tc, t and B. 
The data are presented with statistieal errors (inner bars) and total errors (outer bars). The 
measurements are compared with fits based on a NLO QCD calculation including resummation 
(NLL) and supplemented by power corrections (PC). The fit results are shown as solid lines and 
are extended as dashed lines to those data points which are not included in the QCD fit (see 
Section 5.3). 
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Figure 2: Normalised event shape distributions eorreeted to the hadron level for po and the 
C-parameter. The data are presented with statistieal errors (inner bars) and total errors (outer 
bars). The measurements are eompared with fits based on a NLO QCD ealeulation including 
resummation (NLL) and supplemented by power corrections (PC). The fit results are shown as 
solid lines and are extended as dashed lines to those data points which are not included in the 
QCD fit (see Section 5.3). The symbols and scale factors are defined in Fig. 1. 


determined here using DISASTER++ [29] together with DISPATCH [10]. The parton density 
functions of the proton are taken from CTEQ5M1 [22]. 

The region of F <C 1 is dominated by events with soft and/or collinear parton emissions, 
leading to large perturbative coefficients [30] 


, log^"-^ F 

fZo F 


( 8 ) 


In this region of low F where the bulk of the data lies, the large logarithms need to be re¬ 
summed to all orders in as- Eor the event shape variables considered the resummation has been 
performed by Dasgupta and Salam [10,31] and is available in the DISRESUM package. 

To obtain a good description of the data over the full range of F, it is necessary to add the 
fixed order and the resummed calculations, and to subtract any terms which are counted twice, 
namely the 0(as) and 0{al) terms of the resummed result. There are several valid match¬ 
ing schemes available. Comparing the results from the different schemes leads to a residual 
ambiguity, which is considered in the uncertainty on the prediction. Eor the central values of 
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Figure 3: Mean values of event shape variables eorreeted to the hadron level as a funetion of 
the seale Q. The data, presented with statistieal errors (inner bars) and total errors (outer bars), 
are compared with the results of NLO QCD fits including power corrections (PC). The dashed 
curves show the NLO QCD contribution to the fits. 
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the present fits the modified logR matehing seheme is chosen, with the parameter ppT set to 
two [ 10 ]. 


All event shape variables are subject to non-perturbative effects due to hadronisation, for 
which calculations based on partons need to be corrected. The corrections can be determined 
by using fragmentation models as applied in Monte Carlo programs. For many applications 
these models lead to a reasonable description of the data. However, the interface between 
perturbative and non-perturbative processes is not well defined within these models, and there 
are phenomenological parameters to be tuned. 

An alternative approach has been developed [2] based on the observation that non-pertur¬ 
bative corrections are in general suppressed by powers of (1/Q). For the present event shape 
variables the leading corrections are proportional io 1/Q according to [2]. The power law 
behaviour of the corrections can be described by introducing an effective coupling Oefr, which 
is valid for low scales. In the perturbative region the effective coupling has to coincide with 
the renormalised coupling as{Q). This is conventionally achieved by matching Oeff to at 
a scale fij = 2 GeV. This ansatz results in only one single non-perturbative parameter ao = 
d/c, being the first moment of the effective coupling integrated over the low scale 
region up to the matching scale fij. Power corrections for event shape variables in DIS have 
been calculated to one-loop [32] and two-loop [33] accuracy. 

For the differential distributions the power correction results in a shift of the perturbatively 
calculated distribution [3] 

1 da(F) _ 1 daPQCD(F - apV) 
a dF ~ a If ’ 

where Oi? is of order one and can be calculated perturbatively. For the jet broadening a squeezing 
is applied in addition to the shift, which is absorbed in the coefficient ap [31]. 

The power correction term V is assumed to be universal for all event shape variables. It is 
proportional iol/Q and evaluated to be 


Stt Q 




/do 


, Q K 

In — + — + 1 
F/ Po 


al{Q) 


( 10 ) 


where/3o = 11—2nj/3, Ji" = 67/6—7r^/2—5n/-/9, andn/- = 5 is the number of active flavours. 
The so-called Milan factor AT ~ 1.05 ensures the universality at the two-loop level [33]. 


The simple shift in Eq. 9 cannot be valid over the whole spectrum. At low values of F it 
may be applied only for F 3> apV ~ Pi/Q [3]. Moreover, at large values of F higher order 
corrections are substantial and the NLO calculation is not reliable. 


The mean value of an event shape variable is modified through non-perturbative effects by 
an additive constant [ 1 ] 

{F) = {FY^^^+ aFV, ( 11 ) 

with the same coefficient ap and the same function V as for the distribution (Eq. 10). 

The theory predicts a universal value of ao(F/) of about 0.5 [1]. Eor mean event shape 
values this prediction has been confirmed within 20% in DIS [4] as well as in e+e“ annihila¬ 
tion [7]. Similar conclusions were drawn in an analysis of differential distributions in 
annihilation [7]. 
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5.2 Fit Procedure 


Fits for as and ao are performed via a minimisation using MINUIT [34], whieh for one event 
shape variable is defined as 




( 12 ) 


with V the eovarianee matrix, i and j extending over the bins ineluded in the fit, m the measured 
data points and t(as,aQ) the theory predietion, whieh depends on the free parameters. The 
eovarianee matrix V consists of the sum of the individual covariance matrices from the different 
sources presented in section 3.4: 



(13) 


where l^tat is determined by the unfolding procedure. Only the unfolding and part of the model 
uncertainties are treated as uncorrelated, i.e. these matrices are diagonal. Summing the covari¬ 
ance matrices in this way corresponds to adding the errors in quadrature. The correlated part of 
the model uncertainty is determined as half of the difference in the fitted parameters when the 
fit is repeated, using the program DJANGOH instead of RAPGAP for unfolding. 

5.3 Fits to Spectra 

In order to fit the calculations described in section 5.1 to the measured distributions, the range 
of the distributions to be used has to be specified. The upper bounds in F are given in [10] and 
are motivated by properties of the perturbative calculations. The lower bounds in F are set by 
the behaviour of the power correction, which limits the reliability of the prediction at low Q, 
see Eqs. 9-10. The present analysis makes an effort to extend the fit interval to values of F as 
low as possible. 

For a given event shape, the is first calculated using the highest F and Q bins only. Bins 
corresponding to lower F values and lower scales are successively included provided that the 
of the fit does not increase by more than four units for each additional bin. 

The data are well described by the resummed pQCD calculation supplemented by power 
corrections as shown in Figs. 1 and 2. In many cases even the lowest F bin could be included 
in the fit, which is reasonable since, although the distribution diverges as F —^0, the integral 
within the measured low F bin remains finite. At low Q values the agreement between the 
measurements and the calculation degrades. In this domain the hadronisation effects become 
more important and the simple shift of Eq. 9 is not expected to hold. 

The results of the combined fit for ao and as{mz) are displayed in Eig. 4 and summarised 
in Table 2. The quality of the fits, expressed in terms of per degree of freedom, is found to be 
reasonable. Eor all event shape variables, consistent values for as{mz) and ao are found, with 
a maximum difference of about two standard deviations between r and C. A strong negative 
correlation between as{mz) and ao is observed for all variables. The values of the strong 
coupling as{mz) are in good agreement with the world average [35], shown for comparison as 
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Figure 4: Fit results to the differential distributions of r, B, po, tc and the C'-parameter in 
the (as,ao) plane. The la eontours eorrespond to = Xmin + 1’ ineluding statistieal and 
experimental systematie uneertainties. The value of as (vertieal line) and its uneertainty (shaded 
band) are taken from [35]. 

the shaded band. The non-perturbative parameter oq — 0.5 is confirmed to be universal within 

10 %. 

The theoretical uncertainties on the fitted values of Oq and as{mz) are determined from the 
changes to the results under variation in the procedure as follows: 

• bins with lower boundaries at F = 0 are omitted; 

• the renormalisation scale pr is varied from Q/2 to 2Q; 

• the infrared matching scale pi is varied from 1.5 GeV to 2.5 GeV; 

• the CTEQ proton pdfs are replaced by three versions of the MRST2001 set [36], which 
differ in as{mz) from 0.117 to 0.121; 

• instead of log R the modified M and modified matching schemes [10] are used. 

The fit procedure is repeated for each of these variations separately. 

All fit results, including the individual contributions to the total error, are given in numerical 
form in Table 2. The theoretical error is the dominant contribution to the total uncertainty and 
arises mainly due to the renormalisation scale uncertainty. 
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strong coupling constant a si 

mz) 

event shape variable 

Tc 

T 

B 

Po 

C 

central value 

0.1171 

0.1202 

0.1196 

0.1174 

0.1156 

uncertainties: 






total 

+ 0.0068 

- 0.0062 

± 0.0072 

- 0.0058 

+ 0.0072 

- 0.0064 

± 0.0070 

- 0.0056 

+ 0.0073 

- 0.0054 

total experimental 

± 0.0035 

± 0.0021 

± 0.0014 

± 0.0021 

± 0.0021 

statistical experimental 

± 0.0014 

± 0.0006 

± 0.0004 

± 0.0010 

± 0.0009 

systematic experimental 

± 0.0033 

± 0.0020 

± 0.0013 

± 0.0019 

± 0.0019 

total theoretical 

± 0.0058 

- 0.0051 

± 0.0068 

- 0.0054 

± 0.0071 

- 0.0063 

± 0.0067 

- 0.0052 

± 0.0069 

- 0.0049 

Pr dependence 

± 0.0054 

- 0.0048 

± 0.0058 

- 0.0043 

± 0.0056 

- 0.0044 

± 0.0064 

- 0.0050 

± 0.0069 

- 0.0048 

Pi dependence 

± 0.0002 

- 0.0002 

< 10 -^ 

< 10 -^ 

± 0.0002 

- 0.0002 

< 10 -^ 

fit interval 

± 0.0015 

- 0.0018 

± 0.0007 

- 0.0022 

± 0.0001 

- 0.0009 

± 0.0010 

± 0.0007 

± 0.0003 

- 0.0004 

parton density functions 

± 0.0002 

- 0.0001 

± 0.0003 

- 0.0010 

± 0.0006 

- 0.0007 

± 0.0001 

- 0.0002 

± 0.0002 

- 0.0001 

matching scheme 

± 0.0015 

± 0.0005 

± 0.0036 

± 0.0022 

± 0.0043 

± 0.0043 

± 0.0018 

± 0.0009 

- 0.0005 

- 0.0009 


non perturbative coupling olq^pi = 2 GeV) 

event shape variable 

Tc 

r 

B 

Po 

C 

central value 

0.488 

0.513 

0.519 

0.486 

0.481 

uncertainties: 






total 

± 0.037 

- 0.035 

± 0.034 

- 0.039 

± 0.059 

- 0.049 

± 0.023 

- 0.035 

± 0.028 

- 0.042 

total experimental 

± 0.021 

± 0.025 

± 0.039 

± 0.014 

± 0.008 

statistical experimental 

± 0.009 

± 0.009 

± 0.006 

± 0.006 

± 0.005 

systematic experimental 

± 0.019 

± 0.023 

± 0.038 

± 0.013 

± 0.007 

total theoretical 

± 0.030 

- 0.027 

± 0.022 

- 0.029 

± 0.044 

- 0.029 

± 0.019 

- 0.032 

± 0.026 

- 0.041 

Pr dependence 

± 0.020 

- 0.026 

± 0.018 

- 0.027 

± 0.030 

- 0.028 

± 0.017 

- 0.027 

± 0.022 

- 0.038 

fit interval 

± 0.022 

- 0.007 

± 0.008 

- 0.005 

± 0.030 

± 0.006 

- 0.003 

- 0.016 

± 0.006 

- 0.003 

parton density functions 

± 0.001 

- 0.001 

± 0.006 

± 0.004 

± 0.011 

± 0.003 

± 0.001 

- 0.001 

± 0.001 

- 0.002 

matching scheme 

- 0.005 

- 0.012 

- 0.009 

- 0.023 

± 0.006 

- 0.010 

- 0.006 

- 0.009 

- 0.014 

- 0.014 


correlation coefficient as,ao 

-0.85 

-0.76 

-0.75 

-0.78 

-0.51 

/ d.o.f. (experimental errors) 

1.13 

0.51 

0.81 

1.40 

1.20 


Table 2: Results of simultaneous fits of as{mz) and = 2 GeV) to the distributions of the 
event shape variables tq , t , B, po and C. The statistical and experimental systematic errors as 
well as the theoretical uncertainties are given. 


The good agreement of the results for all event shape variables allows a common set of 
values of and ao to be derived by applying an averaging procedure to the results from 

the individual event shape variables. In this procedure, the minimisation takes into account 
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all experimental and theoretieal errors and the eorrelations between asijnz) and oq as given in 
Table 2. In addition the eorrelations among the observables are eonsidered. The precision on 
the very large correlation coefficients within the group [ t , B} and the group { tq , Po , C*} is not 
sufficient to allow the correlation matrix to be used directly in a minimisation. Instead the 
averaging is performed in two steps. Firstly, the results within these two groups are combined, 
neglecting the correlations among the observables. The smallest uncertainty of the contributing 
measurements is conservatively taken as the uncertainty of the group average. Secondly, the 
two group averages, being only moderately correlated, are combined using the corresponding 
correlation matrix. 


The averaging procedure results in: 

as{mz) = 0.1198 ± 0.0013 (exp) ^o;oo 43 (theo) , 
ao = 0.476 ± 0.008 (exp) lo;o 59 (theo) , 


(14) 


with a fit quality of x^/d.o.f. = 4.9/2. Here the theoretical error is derived from the renormali¬ 
sation scale uncertainty. Note that the combined value of is lower than the individual values, 
due to the negative correlations between ao and as{mz)- 

If instead the correlations between the two groups of event shape observables are neglected, 
a consistent result is obtained. 


5.4 Running of Strong Coupling cXs{Q) 

In the previous section it was shown that the concept of power corrections provides a good 
description of hadronisation effects in the differential event shape distributions. Alternatively, 
one may assume the validity of the power corrections and investigate the scale dependence of 
the strong coupling as{Q). For each event shape variable a Q independent oq parameter and an 
as{Q) for each Q bin are fitted. 

The fitted values of as{Q) are presented in Fig. 5. The running of the strong coupling is 
clearly observed for each of the event shape variables over a wide scale range between {Q) = 
15 GeV and {Q) = 116 GeV. The numerical values with experimental errors are given in 
Table 3. The theoretical uncertainties are of similar size to those given in Table 2. 

The fit results of the different event shapes are compatible with each other and again may 
be combined. The correlations among the observables are taken into account in a two-step 
procedure as described in the previous section. The averaged as{Q) values are displayed as a 
function of Q in Fig. 6 and listed in Table 3. A fit of the renormalisation group equation to the 
measured as{Q) yields 

Oisimz) = 0.1178 ± 0.0015 (exp) lQ;QQg} (theo), (15) 

with x^/ndf = 8.3/6. This value is in good agreement with the two-parameter fit result quoted 
in Eq. 14, though the scale uncertainty is somewhat larger here. The difference with respect to 
the result quoted in Eq. 14 is that here individual power correction parameters oq are associated 
to each observable. 
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Q/ 







strong 

coupling 

constant 

OLsiQ) 







GeV 

Tc 

T 

B 

Po 

C 

(a 

m) 

15 

0.155 

± 

0.007 

0.171 

± 

0.008 

0.164 

± 

0.004 

0.162 

± 

0.006 

0.162 

± 

0.006 

0.163 

± 

0.004 

18 

0.149 

± 

0.007 

0.165 

± 

0.005 

0.160 

± 

0.003 

0.158 

± 

0.005 

0.157 

± 

0.006 

0.159 

± 

0.003 

24 

0.150 

± 

0.005 

0.148 

± 

0.006 

0.151 

± 

0.003 

0.147 

± 

0.004 

0.147 

± 

0.004 

0.148 

± 

0.003 

37 

0.133 

± 

0.006 

0.139 

± 

0.004 

0.141 

± 

0.003 

0.137 

± 

0.004 

0.130 

± 

0.004 

0.136 

± 

0.003 

58 

0.128 

± 

0.006 

0.127 

± 

0.005 

0.129 

± 

0.005 

0.127 

± 

0.006 

0.131 

± 

0.006 

0.128 

± 

0.005 

81 

0.110 

± 

0.005 

0.107 

± 

0.009 

0.108 

± 

0.006 

0.110 

± 

0.006 

0.104 

± 

0.006 

0.108 

± 

0.004 

116 

0.090 

± 

0.018 

0.138 

± 

0.047 

0.123 

± 

0.027 

0.112 

± 

0.025 

0.122 

± 

0.011 

0.114 

± 

0.011 

ao 

0.496 

± 

0.014 

0.515 

± 

0.020 

0.516 

± 

0.034 

0.484 

± 

0.011 

0.477 

± 

0.007 



Table 3: Fitted values of as{Q) as determined from the differential distributions of the event 
shape variables tq, r, B, po and C and the averaged values {as{Q)) ; the fitted parameters cto 
are given in the last line. 


5.5 Fits to Mean Values 

The mean values of the event shape variables, presented as a funetion of the scale Q in Fig. 3, 
are also subjected to QCD fits. For this application the resummed calculation can not be used, 
because of difficulties at high values of the event shape variable, which are related to sub-leading 
logarithms and the matching procedure [28]. These regions are excluded from the fits to the full 
spectra, but by definition contribute to the mean values. Hence the theoretical prediction for the 
mean values is solely based on a NLO calculation, supplemented by power corrections. 

The results of the fits are displayed in Fig. 3 and the fitted values of as{rnz) and Oq are 
shown in Fig. 7. The non-perturbative parameters cluster around a common value of ao ~ 
0.45 — 0.50, thus supporting universality at the 10% level. The fitted values for as{mz) exhibit 
a rather large spread. 

Compared to the results obtained from the event shape distributions shown in Fig. 4, higher 
values of are found using tc, Po and in particular the C-parameter, and a somewhat lower 
value is obtained from B. The analyses of both the event shape distributions and mean values are 
based on the same data and there is no experimental reason why they should lead to different 
results. Therefore, it is likely that the different theoretical treatment for mean values is the 
source of the observed deviations. 

A dedicated study of mean values has been published previously by the HI collaboration [4] 
for 7 < Q < 100 GeV. As mentioned in section 4 the current and previous measurements are in 
agreement in the phase space of overlap. The results of both QCD analyses are also consistent 
with each other. The previous analysis shows a larger sensitivity to the parameters as and ao 
because it includes data at lower Q. This leads to substantially reduced error ellipses compared 
to those of Fig. 7, in particular for the (r) observable. 

Differences between sets of QCD parameters a^ and ao determined in differential distribu¬ 
tions and mean values, similar to those reported here, have been observed in annihilation 
by DELPHI [7]. Discrepancies have also been found by other experiments [6, 8]. In general, 
similar values of ao — 0.5 are found for event shape variables in e+e“ scattering, but with 
a much larger spread than in deep-inelastic scattering, the jet broadening results in particular 
being different from the others. 
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Figure 5: The strong coupling as a function of the scale Q. The individual fit results, shown 
as points with error bars, are obtained from fits to the differential distributions in re, r, 5, po and 
C within each Q bin. The errors represent the total experimental uncertainties. For each event 
shape observable a value of as{mz) is indicated in the plot, determined from a fit to the as{Q) 
results using the QCD renormalisation group equation. The corresponding fit curves are shown 
as full lines. The shaded bands represent the uncertainties on as{Q) from renormalisation scale 
variations. 

6 Conclusions 

Accurate measurements of event shape variables in deep-inelastic ep scattering are presented 
based on 106 pb“^ of data with four-momentum transfer Q ranging between 14 GeV and 
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Figure 6: The strong eoupling Og as a function of the scale Q from an average of the results 
obtained by fitting the differential event shape distributions. The errors represent the total ex¬ 
perimental uncertainties. A value of as{mz) is indicated in the plot, determined from a fit to 
the as{Q) results using the QCD renormalisation group equation. The fit curve is shown as 
the full line. The inner (outer) shaded band represents the uncertainty of the fitted as{Q) from 
experimental errors (the renormalisation scale variation). 


200 GeV. Resummed perturbative QCD predictions together with power corrections give good 
descriptions of the spectra of the observables thrust, jet broadening, jet mass and G-parameter. 
The use of resummed calculations extends the good description to low values of the event shape 
variables, corresponding to pencil like configurations. 

The results of a two-parameter fit of the strong coupling constant ag and the effective non- 
perturbative coupling ao for the various event shape observables are consistent with each other. 
The values for a* agree with the world average. The parameter ao, which accounts for hadro- 
nisation, is consistently found to be 0.5 within 10%, in good agreement with theoretical expec¬ 
tation. A combined analysis of all event shape variables yields 

o^simz) = 0.1198 ± 0.0013 (exp) lo;oo 43 (theo) , 
ao = 0.476 ± 0.008 (exp) lo;o 59 (theo) , 

where the theoretical error is derived from the renormalisation scale uncertainty. Relaxing the 
requirement of a common value of ao, the data are used to investigate the scale dependence of 
the strong coupling over a wide range ofQ = 15 — 116 GeV. The running of ag{Q) is clearly 
observed for each event shape variable, in accordance with the expected evolution. Combining 
the results of all variables leads to ag{mz) = 0.1178 ± 0.0015 (exp) lo;ooeJ (theo). The errors 
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Figure 7: Results of fits to the mean values of r, B, pQ, tq and the C-parameter in the (a^, ao) 
plane. The la eontours eorrespond to = Xmin + 1’ ineluding statistieal and experimental 
systematie uneertainties. The value of as (vertieal line) and its uneertainty (shaded band) are 
taken from [35]. Note the enlarged seale eompared to Fig. 4. 

are dominated by the renormalisation scale uncertainty, which suggests that missing higher 
order terms in the perturbative calculation are important. 

In the analysis of the event shape means the results for as{mz) and oq are less accurate 
than those obtained in the fits to the distributions. While the non-perturbative parameters a^ 
cluster again around a common value of 0.5, the results for the strong coupling ag obtained 
from the five event shape means exhibit a spread considerably larger than is expected from 
their individual uncertainties. These discrepancies may be related to an insufficient theoretical 
treatment, which lacks resummed calculations. 

The observed universality of ao for both distributions and mean values of the event shape 
variables supports the concept of power corrections. Hence, it can be considered as an appro¬ 
priate alternative to conventional models for the description of hadronisation effects for event 
shape variables. 
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